Pseudorandom noise code-based technique for thin cloud discrimination with C0 2 and 0 2 

absorption measurements 


Joel F. Campbell, Narasimha S. Prasad, Michael A. Flood 

NASA Langley Research Center, 5 N. Dryden St., Hampton, VA 23681 

joel.f.campbell@nasa.gov 


Abstract 

NASA Langley Research Center is working on a continuous wave (CW) laser based remote sensing 
scheme for the detection of C0 2 and 0 2 from space based platforms suitable for ACTIVE SENSING 
OF C02 EMISSIONS OVER NIGHTS, DAYS, AND SEASONS (ASCENDS) mission. ASCENDS 
is a future space-based mission to determine the global distribution of sources and si nk s of 
atmospheric carbon dioxide (C0 2 ). A unique, multi-frequency, intensity modulated CW (IMCW) 
laser absorption spectrometer (LAS) operating at 1.57 micron for C0 2 sensing has been developed. 
Effective aerosol and cloud discrimination techniques are being investigated in order to determine 
concentration values with accuracies less than 0.3%. In this paper, we discuss the demonstration of a 
pseudo noise (PN) code based technique for cloud and aerosol discrimination applications. The 
possibility of using maximum length (ML)-sequences for range and absorption measurements is 
investigated. A simple model for accomplishing this objective is formulated, Proof-of-concept 
experiments carried out using SONAR based LIDAR simulator that was built using simple audio 
hardware provided promising results for extension into optical wavelengths. 
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INTRODUCTION 

The National Research Council’s (NRC) Decadal Survey (DS) of Earth Science and Applications 
from Space has identified the Active Sensing of C0 2 Emissions over Nights, Days, and Seasons 
(ASCENDS) as an important Tier II space-based atmospheric science mission. The C0 2 mixing 
ratio needs to be measured to a precision of 0.5 percent of background or better (slightly less than 2 
ppm) at 100-km horizontal resolution overland and 200-km resolution over oceans. To meet this 
goal, the ASCENDS mission requires simultaneous laser remote sensing of C0 2 and 0 2 in order to 
convert C0 2 column number densities to average column C0 2 mixing ratios (XC0 2 ). As such, the 
C0 2 column number density and the 0 2 column number density will be utilized to derive the average 
XC0 2 column. The anticipated benefits of ASCENDS as discussed in the decadal survey include (a) 
quantification of global spatial distribution of atmospheric C0 2 on scales of weather models as well 
as terrestrial and oceanic sources and sinks of C0 2 during day/night over all seasons, and 
(b)establish a scientific basis for future projections of C0 2 sources and si nk s through data-driven 
enhancements of Earth-system process modeling. 



Accordingly, NASA Langley Research Center (LaRC), with its partners, is working on a CW laser 
absorption spectrometer based remote sensing scheme operating in the 1.57 micron spectral band for 
the detection of CCL and 1.26 micron spectral band for the detection of CL For concentration 
detennination, differential absorption lidar (DIAL) scheme at a selected transition for each gas is 
adopted. Hence, two are more wavelength at a known transition are utilized. The 1.5 micron spectral 
band lying in the telecom region was chosen for CCL detection due to spectroscopic properties 
combined with architectural advantages. The selection of 1.26 micron band for CL detection that is 
close to 1.57 micron spectral band instead of 0.76 micron spectral band is anticipated to provide 
processing advantages. Within these spectral bands, an optimal transition that is least sensitive to 
environmental parameters for each of these trace gases is selected. With limitations on CW laser 
powers, the remote sensing systems may have to operate in the photon starved regions. However, an 
Intensity modulated continuous wave (IMCW) technique is adopted to achieve highest possible SNR 
for processing return signals to minimize the influence of various noise sources within the filter 
bandwidths for each channel. 

Unlike pulsed laser based DIAL systems, CW laser based remote sensing systems are affected 
by cloud and aerosol interferences. The reason for this is because the optical depth, which is a 
measure opaqueness of the atmosphere at a particular optical path length, is a function of the 
distance to the target in addition to the concentration of whatever chemical species one is trying to 
measure at the wavelength of the laser. One could use an additional laser altimeter to measure 
distance and that has been done in the past with some degree of success. However, in situations 
where there are thin clouds the distance one uses in the optical depth calculation is ambiguous 
because we have a return signal from the ground and thin cloud at the same time. The technique 
presented here gives one the ability to discriminate between the return from a cloud or the ground in 
a natural way by calculating a return profile as a function of distance in a manner analogous to pulse 
lidar. In essence, one uses a modulation scheme in combination with signal processing to convert 
that modulation pattern into a pulse. One laser wavelength is tuned to an absorption line of CCL 
which is the on-line wavelength, for instance, while the second laser wavelength is tuned to a 
wavelength just off that wavelength, which we call the off-line wavelength. One separates the two 
wavelengths by modulating each wavelength with an orthogonal sequence because a photo detector 
is not able to distinguish between such close wavelengths. Signal processing is used to separate out 
each wavelength and obtain a pulse profile for each wavelength. By taking the ratio of the pulse 
amplitudes for each wavelength we may calculate the CO 2 concentration from the optical depth and 
the distance. Accordingly, various realizable aerosol and cloud discrimination techniques are being 
investigated for determining required accuracies less than 0.5%. In this paper, the adaptation of a 
Pseudorandom Noise (PN) code based technique to discriminate returns from clouds, aerosols and 
ground is discussed 1 . This is followed by proof-of-concept experiments carried out using acoustic 
frequencies. 

Resolution and Range 

The minimum resolution is given by 
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where c is the speed of light and B R is the bit rate of the PN code. The maximum range is given by 

R = Nr (2) 

where N is the code length. For example, for a PN code of length 255 bits and a bit rate of 50 kHz, 
the maximum range would be 764 km, which is the range we would need for a satellite. However, 
there are advantages to making the code longer as discussed in Section below. 

PN CODES 


PN codes are well kn own in the communications devices for their use in encoding communication 
channels such as in spread spectrum communications. One of the simplest PN codes is the m- 
sequence. Furthermore, maximum length ml-sequences are basis of many other PN codes. We use 
(ml) sequences. Besides radar and sonar, the ml-sequences have been used in Lidar applications for 
measuring range ’ . Pulse lidars inherently provide range information. However, in the case of CW 
laser based remote sensing systems, range information on clouds and hard targets can be obtained by 
transmitting amplitude modulated PN code. Range profiles can be extracted by cross correlation 
between the reflected signal and the transmitted code. The correlation peaks of range profiles will 
indicate the location of scattering centers such as clouds and ground. 

New codes that are orthogonal to each other over a limited operational range can be generated. 
This gives us the ability to design a multichannel DIAL systems. As such, it allows one to measure 
relative absorptions of two or more laser wavelengths at a particular range to detect specific 
atmospheric species. Besides m-sequences 4 ’ 5 other codes are also useful for ranging 6 . However, the 
use of time shifted PN codes to modulate multiple laser wavelengths for measuring absorption at 
specific ranges in order to discriminate the return of a cloud from the ground is a very effective and 
simple solution. For our application, the m-sequences are especially advantageous because they have 
very good autocorrelation properties. 

The ML sequences can be represented in a number of different ways. One popular method is by 
using linear feedback registers with modulo 2 additions 7 which is shown in Figure 1. 



Figure 1 . Shift register representation for m-sequence codes. 
Another representation is using a generator polynomial 


G = g m x m + + g m _ 2 * m - 2 + •••+ SA 2 + &* 1 + So 
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Here the g’s can be 1 or 0 and the sums are done using modulo two addition. However, g m =g 0 =l. 


These sequences are 2 m - 1 in length and only specific polynomials can be used for a particular length. 
For m=8 which corresponds to a code length of 2 8 -l=255, an allowed polynomial is, 

G = + x 6 + x 5 + x 1 +\ (4) 

In order to generate the individual bits, we can represent this by the recursion relation 

X j+& = Xj ®X j+2 ®X j+5 ®X j+6 , (5) 

where X is the maximum length sequence and j is the index. For the initial values of X(, - X-/ we may 
choose any sequence of Os and Is we wish as the seed as long as they are not all 0 or all 1. For the 
seed (1,0, 1,0, 1,1, 1,1) we generate the repeating sequence 


(1, 0, 1, 0, 1, 1, 1, 1, 0, 0, 1, 1, 1, 0, 1, 0, 0, 0, 0, 1, 0, 1, 0, 1, 1, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 
0 , 0 , 1 , 1 , 0 , 0 , 1 , 0 , 0 , 0 , 1 , 1 , 0 , 0 , 0 , 0 , 1 , 1 , 0 , 1 , 1 , 1 , 1 , 1 , 1 , 0 , 1 , 1 , 1 , 0 , 0 , 0 , 0 , 1 , 0 , 0 , 0 , 0 , 0 , 1 , 

0 , 0 , 1 , 0 , 1 , 0 , 1 , 0 , 0 , 1 , 0 , 1 , 1 , 1 , 1 , 1 , 0 , 0 , 0 , 0 , 0 , 0 , 1 , 1 , 1 , 0 , 0 , 1 , 1 , 0 , 0 , 0 , 1 , 1 , 0 , 1 , 0 , 1 , 0 , 0 , 

0 , 0 , 0 , 0 , 0 , 1, 0 , 1, 1, 1, 0 , 1, 1, 1, 1, 0 , 1, 1, 0 , 0 , 1, 1, 1, 1, 1, 1, 1, 1, 0 , 0 , 1, 0 , 1, 1, 0 , 1, 0 , 1, 1, 0 , 

1 , 1 , 1 , 0 , 1 , 0 , 1 , 0 , 1 , 0 , 1 , 1 , 1 , 0 , 0 , 1 , 0 , 0 , 1 , 1 , 0 , 1 , 1 , 0 , 1 , 0 , 0 , 1 , 1 , 0 , 0 , 1 , 1 , 0 , 1 , 0 , 0 , 0 , 1 , 1 , 

1 , 0 , 1 , 1 , 0 , 1 , 1 , 0 , 0 , 0 , 1 , 0 , 0 , 0 , 1 , 0 , 0 , 1 , 1 , 1 , 1 , 0 , 1 , 0 , 0 , 1 , 0 , 0 , 1 , 0 , 0 , 0 , 0 , 1 , 1 , 1 , 1 , 0 , 0 , 0 , 

1 , 0 , 1 , 0 , 0 , 1 , 1 , 1 , 0 , 0 , 0 , 1 , 1 , 1 , 1 ) 


Correlation properties 

As mentioned earlier, m-sequence codes have very good autocorrelation properties. The cross 
correlation between itself and a shifted version of itself is N (code length) when they are in sync and 
-1 if they are not. It is also understood that the codes are AC coupled before any correlation is done 
by changing all O’ s in the code by -1. There are a number of ways of computing the correlation. One 
is by computing a table of 

N 

R[A,B] = Y, A[m]B[m + n] (6) 

m = 1 


Note that some people define this differently by dividing by the length. This works but it requires 
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something of the order of N“ steps. A better and faster way is to use Fourier transforms. In this case, 

R[A,B] = F- l [F[Aj F[B]] (7) 

where the * indicates the complex conjugate and F is the FFT. This is superior because it only takes 
something of the order of N log N steps. As an example we take the case in the previous section and 



compute the cross correlation between it and a code shifted 128 places to the right and the result is 
shown in Figure 2. 



Figure 2. Cross correlation of code and shifted code computed with FFT. 


Multiple channels 

In order for this to be useful for atmospheric species detection it is advantageous to be able to 
separate out multiple transmitted wavelengths. Each wavelength is modulated with a different 
orthogonal code. We may take advantage of two facts to use m-sequences (for instance) to construct 
a system with two or more channels, where each channel represents a different wavelength that is 
modulated by a different code. Since m-sequences have very good autocorrelation properties shifted 
versions of them will cross correlate well except for a small DC component where they are 
uncorrelated and the big spike where they correlate. The second thing we have with this application 
is when the LIDAR is looking down from space there will be a final hard target (the ground) below 
which there will be no further returns possible. That opens up some interesting possibilities. If we 
know there will be no further returns past a certain distance we can size a system with twice the 
range we need and use that further range to add a second shifted version of that same code that is 
shifted by half the doubled range. To get twice the range we simply use a PN code twice as long 
while maintaining the same code bit rate. As long as there are no returns beyond that half range there 
will be no interference between the channels. The first half of the cross-correlation gives the range 
info for the first unshifted channel and the second half of the autocorrelation gives the range info for 
the shifted second channel and we can do that in a single step. By doubling the length of the code we 
also increase the SNR as we see from the previous section. We can make the code as long as we 
need to add more channels. 

In the sample below we use a code length of 5 1 1 with a sample rate of 500000 samples/sec and a 
code bit rate of 50000 bits/sec. This gives us an over-sampling of M=10 and gives us a total range of 
51 1*3*1 0 A 5/(2* 50000)= 1530 km. To generate the code we use a code with the generating function 
G = x 9 + x 1 + x 3 + x 1 + 1 , a seed of (1,0, 1,0, 1,1, 1,1,1) and a shifted version of that code shifted 256 
places to the right. The result of that and lOx noise is shown in Figure 3. 
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Figure 3. Example of using two codes in a double ranged system with two channels. 

Signal-to-noise ratio estimates 

We first consider white noise. It is possible to estimate the signal to noise ratio improvement before 
and after correlation. Let Sr be the sample rate and Br be the PN code bit rate. Then number of 
samples per code bit is 


M = — . (8) 

B r 

Oversampling the code by a factor of M increases the length of the code by the same factor so the 
code length becomes MN. 

Lets suppose we have an AC coupled signal which we represent by 

S = r) + a pn A (9) 

where 77 is the noise, a is the attenuation factor and 


pnf = 2PN i+A - 1 (10) 

The term PN is the PN code represented by l’s and 0’s and A represents the phase shift and pn 
contains l's and - l’s. The standard deviation of the noise we represent by 



- n ave f 
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where we have summed over the total length of the code. After perfonning the correlation 
calculation, we see that the signal strength increases by a factor of MN where they are correlated. 
We now ask the question what correlation does to the standard deviation of the noise. The 
correlation function of the total signal is, 

MN MN MN 

R U) = X('f + a p n t )p n lj p r lj + ®Lp n t p n lj = Vj ,+ Sj 1 (12) 

i = 1 i=l i = 1 

The question now is to find the standard deviation of the transformed noise 7] ’ which we assume to 
have 0 mean (AC coupled). We first find the variance by 

i MN 

(<A=(('if)=] j ^S(v) < 13 > 


To find the sum we write 

MN MN MN 

(Vj ’)" = Efim P n l + ^n P lt n +J = ^ V n P>N + j P'Nj ( 14 ) 

m= 1 n = 1 n,m 


so that 
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j = 1 n,m 


If the noise has 0 mean, the only part of the sum that will contribute is when n=m so 

i MN MN i MN MN i MN 

(<7') 2 = — HvilpnlJ = — 1 1 r£ =— t MNo 2 -MNo 2 
V ’ MN%^ ' n[ m+,) MN MN 

In the above since pn is either a 1 or -1, its square is always 1. The result is that 

o' = \JmNg . 


(16) 


(17) 


The new correlated snr in terms of the original snr is 
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( 18 ) 


Expressed differently in tenns of previously defined variables we find 
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It should be stressed that this was derived assuming noise with 0 mean. If the DC offset is not 
subtracted from the signal before processing we would get a very different result because the PN 
code has a small DC offset for the simple reason it has an odd length. Alternatively one may subtract 
the offset from the code. The correlation properties would be similar and the resulting code would be 
orthogonal to any DC bias in the data. 

Photo detector SNR 

If the system is lidar we must use the SNR equation for photodiodes. As an example for two channel 
PN code modulation we use an on-line and a time shifted off-line code to modulate the on-line and 
off-line wavelengths. The SNR for such a signal is. 


SNR = 7 l -tL = . (20) 

VMCWO W 

Here Is represents the signal current and In represents the equivalent noise current from all noise 
sources, whereas i is the AC coupled part of I. The variables q and A/ are the electron charge and 
noise bandwidth. For the received PN signal the SNR for each channel would be 
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( 21 ) 


Here (P sl ) and (P s2 ) are the average on and off-line signal power, whereas (P v ) and 01 are the 

equivalent total noise power and responsivity. We use the RMS power in the numerator because 
only the AC coupled part of the signal is used in the processing. The variable p represents the AC 
coupled part of P. We first write 

pn = 2PN-\. (22) 

Since (PN) = ( N + 1) / (2 N ) , we find 


P sn = ( P s,,)j^- l ( l + P n )- 


(23) 


The AC coupled part of P is 
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The 1/N term is the slight dc offset of pn. From this we find 
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Therefore, 
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We may now use the results of Equation 18 to find the SNR after correlation. If we subtract the bias 
from the noise before correlation we find the new SNR after correlation is 
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where n and «' take on the values of 1 or 2 and n + n' . The extra MN term under the square root 
comes about because that part of the noise is correlated with the signal. This compares favorably 
with what has been reported elsewhere" if the differences in peak versus average power, RMS 
power, over sampling, and number of channels are taken into account. 
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Figure 4. Pulse lidar system for on-line/off-line measurement 


Comparison with pulse 


For a pulse lidar system as shown in Figure 4 the SNR is 


SNR pidse = 


SIP, 


pulse Sn 




( 27 ) 


for one complete on-line/off-line cycle. In order to match the resolution for pulse we have 
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If we make the sample times the same for a complete cycle we have N = 2T B r = IT / r . If we 
make the noise bandwidth the same and make the average signal power for the PN system the same 
as the pulse system multiplied by a factor of p then by Equation 26 the PN SNR is 
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For the same noise bandwidth and noise levels we find 
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In the strong signal limit we have 
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In the weak signal limit we find 


SNR pulse F V N + 1 V 2T 


(32) 


The strong signal limit is a case we would not likely see since the detector shot noise will typically 
be less than other noise sources. As a result we are left with the weak signal limit. We may use this 



to optimize the CW system for the best possible performance in terms of sample rate and average 
power level. Of course there are other considerations as well. In any real system one must consider 
the specific implementation when making such a comparison. One obvious advantage of the CW 
system is it can make a simultaneous on-line/off-line measurement whereas with the pulse system 
the measurement is sequential. This can potentially make the pulse system less accurate in 
detennining the optical depth ratios. 


EXPERIMENTAL RESULTS AND DISCUSSION 
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To test these concepts we implement the system shown in Figure 5 and 6 in audio ' . Because an 
audio speaker is not capable of reproducing the very low frequencies required by the PN code, one 
must first put the PN code on a carrier, then modulate that carrier with the PN code. Here the 
resolution and range are represented by Equations 2 and 3 except that c becomes the speed of sound 
in this case. 

o 

In the initial development, this system was first implemented in Mathematica then in Labview. 
In order to accomplish this, special Labview audio drivers written by Christian Zeitnitz 10 were 
utilized. This increased both the performance and functionality of the Labview audio interface. A 
Samson Servo 201a audio amp was used for the audio amplifier. This had a bandwidth of 65 kHz 
and is designed for accurate sound reproduction. The audio card was an M-Audio 1814 with a 1394a 
interface. This card is capable of sample rates of up to 192 kHz. The software display is shown in 
Figure 7. 



Figure 5. Audio implementation 
for single channel audio ranging 
system. All 

modulation/ demodulation, filtering, 
and cross correlation functions 
were implemented in Labview 
software. 


A dual channel implementation is also possible which could potentially be used for differential 
absorption measurements. This is analogous to the proposed LIDAR implementation and can 
potentially be used for differential absorption measurements. 



Figure 6. Dual channel 
implementation of audio ranging 
system can potentially be used for 
differential absorption 
measurements. This is analogous 
to how it could potentially be 
implemented in LIDAR. 





Figure 7. Labview interface and display for the dual channel implementation. 


The proposition is to use two or more PN codes where each is shifted in time in a LIDAR 
implementation. The shifting of the codes is done such that there is equal spacing. For instance with 
two codes of length 511, the second is shifted to the right 256 places with respect to the first code. 
The shifting is done in such a way that as the end of the code is shifted past the boundary it wraps to 
the beginning. The maximum distance to the farthest target shall be 255 in that case. With four codes 
of length 511, each is shifted 128 places to the right of the preceding code. The maximum distance 
to the target in that case shall be 127. In this way one may have multiple channels for a given code 
length. The only restriction is that the LIDAR/RADAR/SONAR must be pointed at targets with a 
hard surface (such as the surface of the earth) no farther that the distance to the first code. Without 
that one will have range-wrapping issues. In addition to this one should either subtract the average 
from the code or from the data prior to processing. This helps eliminate noise/bias issues. 


CONCLUSIONS 

In this paper, a novel method based on PN codes to discriminate clouds and aerosols from ground 
returns suitable for CW laser based remote sensing systems for space based platforms is discussed. 
The intended use of this technique is in the ongoing for ACTIVE SENSING OF C02 EMISSIONS 
OVER NIGHTS, DAYS, AND SEASONS (ASCENDS) program effort at NASA Langley Research 
Center. It is shown that time-shifted PN codes would be valuable for obtaining range information 
of interfering clouds and aerosols so that returns from ground or topographic background can be 
discriminated during data retrieval processes. Proof-of-concept prototype operating at audio 
frequencies has been successfully demonstrated. Plans are underway for demonstration at optical 
frequencies. 
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